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Abstract. We discuss the effects of a trapping space-dependent potential on the 

critical dynamics of lattice gas models. Scaling arguments provide a dynamic trap-size 
^ ' scaling framework to describe how critical dynamics develops in the large trap-size 

^5 . limit. We present numerical results for the relaxational dynamics of a two-dimensional 

lattice gas (Ising) model in the presence of a harmonic trap, which support the dynamic 

trap-size scaling scenario. 
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Statistical systems are generally inhomogeneous in nature, while homogeneous 
systems are often an ideal limit of experimental conditions. Critical phenomena, 
characterized by the development of correlations with diverging length and time scales, 
^ • are usually considered assuming homogeneous systems. However, the emergence of 



O 



critical correlations is also observed in inhomogeneous systems. Particularly interesting 

physical systems are interacting particles constrained within a limited region of space 

by an external potential. This is a common feature of the experimental realizations of 
^ ' 1—1 

ly-^ ■ Bose- Einstein condensations in diluted atomic vapors [1] and of optical lattices of cold 

atoms [2] . Experimental evidence [3] of a critical behavior has been reported for a three- 

O ' dimensional trapped Bose gas, observing an increasing correlation length compatible 

t^^ ■ with a continuous transition, although the system is made inhomogeneous by the 

■ confining force. Experimental evidences of finite-T Kosterlitz-Thouless transitions |4] 

in two-dimensional trapped Bose atomic gases have been provided in [HI E] [TJ El [9]. 

However, the inhomogeneity due to the trapping potential is expected to strongly 

^ ■ affects the phenomenology of continuous phase transitions observed in the absence of 

d • a trap. For example, correlation functions of the critical modes are not expected to 

develop a diverging length scale in a trap. Therefore, a theoretical description of the 

critical correlations in systems subject to confining potentials, and of how they unfold 

approaching the transition point, is of great importance for experimental investigations 

of the critical behavior of systems of trapped interacting particles. 

We consider a trapping potential 

V{r)=vP\r\P={\f^/iy, (1) 

where v and p are positive constants and / = v~^ is the trap size, coupled to the particle 
number. Harmonic potentials, i.e., p = 2, are usually realized in experiments. The effect 
of the trapping potential is to effectively vary the local value of the chemical potential, 
so that the particles cannot run away. Let us consider the case in which the system 
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parameters, such as temperature, pressure and chemical potential, are tuned to values 
corresponding to the critical regime of the unconfined system. The critical behavior 
gets distorted by the trap, although it gives rise to universal effects in the large trap- 
size limit, controlled by the universality class of the phase transition of the unconfined 
system. The corresponding scaling regime can be cast in the form of a trap-size scaling 
(TSS) [TOllTT], resembling the finite-size scaling theory for homogeneous systems, but 
characterized by a further nontrivial trap critical exponent 6, which describes how the 
length scale ^ of the critical modes depends on the trap size at criticality, i.e., ^ ~ l^, 
and which can be estimated using renormalization-group (RG) arguments. 

In this paper we investigate how critical dynamics develop in trapped systems at 
thermal continuous transitions, i.e., c/asszca/ transition driven by thermal fluctuations. 
We extend the TSS Ansatz of [10], to allow for the time dependence of the critical 
modes, which depends on the static universality class, but also on the general features 
of the dynamics, analogously to homogeneous systems [12] . 

In a standard scenario for a continuous transition, see, e.g., [I3], the critical 
behavior of a rf- dimensional system is characterized by two relevant parameters Ut and 
Uh, which may be associated with the temperature T, i.e., Ut ~ T/Tc — 1, and the 
external field h coupled to the order parameter, whose RG dimensions are yt = l/iy and 
yh = {d + 2 — If])/ 2. The presence of a trap of size / generally induces a further length 
scale ^trap ~ l^ , which must be taken into account to describe the critical correlations. 
Within the TSS framework [lOl IT], see also [Tail5l|16l[I71[l8l[l9l[20l[2ll[22], the 
scaling law of the singular part of the free energy density around the center of the trap 
can be written as 

where 9 is the trap exponent. At the critical point ut = 0, the TSS of the correlation 
function of the order parameter (p is 

G{x, y) ^ {<P{x)4>{y)), = V'^'-^^^^ f {xV^ yV'), (3) 

thus implying ^ ~ /^ for its length scale. Finite size effects, due to a finite volume L'^, 
can be taken into account by adding a further dependence on Ll~^ in the above scaling 
Ansatz [19] . 

The universal scaling behavior of the critical dynamics is essentially determined 
by the dynamic universality class, which depends on a few general properties of the 
dynamics, for example, whether there are conserved quantities during the time evolution 
(see [12] for a list of dynamic universality classes) . The scaling behavior of the correlation 
functions requires a further dynamic critical exponent 2;, which provides the power- 
law relation between the critical time scale r and the diverging correlation length 
T ~ ^^ ~ M^^^, thus implying a diverging time scale at the critical point. For example, 
the equilibrium time dependence of the correlation function of the order parameter 
0(x, t) is expected to scale as 

G(xi, X2; ti, t2) = (0(xi, ti)0(x2, t2))c = e~^'~'+"V[(a;2-xi)/e, (t2-ti)/r] (4) 
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Again, the presence of the trap drastically affects the critical dynamics; in particular, 
the time correlations are not expected to develop a diverging time scale. Within the 
TSS framework, the time dependence of the dynamic scaling behavior is expected to 
enter through a further dependence on the scaling variable tl~^^. Thus, at the critical 
point Ut = 0, the scaling behavior at equilibrium of the two-point correlation function 
dl]) is expected to be 

G{x^,X2■, t,M) = ^^('^-2+'')/,[x2/-^ x,l-\ (t2 - ti)/-'1 (5) 

which implies that the corresponding time scale r behaves as 

r ~ f' (6) 

at the critical point. Analogous scaling behaviors have been put forward to describe 
the temperature and time dependence at T = quantum transitions [HI [231 El]- The 
dynamic TSS Ansatz can also be extended to take into account off-equilibrium dynamics. 
In order to check the dynamic TSS scenario, we consider the lattice gas model 
defined by the Hamiltonian 

^Lgas = -4 J ^ p^pj - l^^Pi + ^ '2V{ri)pi (7) 

where the first sum runs over the nearest-neighbor sites of a square lattice, pi = 0, 1 
depending if the site is empty or occupied by the particles, /x is the chemical potential, 
and V{r) is the potential ([I]) which vanishes at the center of the trap. Far from the 
origin the potential diverges, thus the expectation value of the particle number tend to 
vanish, indeed (px) ~ exp[— 2l^(x)] at large distance from the center, and therefore the 
particles are trapped. The lattice gas model (171) can be exactly mapped to a standard 
Ising model, by replacing pj = (1 — Sj)/2, obtaining 

7^ = -J^SiS, -/i^Si-^l^(ri)si (8) 

(ij) i i 

where Sj = ±1 and h = 8J + fi/2. In this picture the external potential plays the role 
of a space-dependent magnetic field. In the absence of the trap, the square-lattice Ising 
model (JH]) shows a critical behavior characterized by the critical exponents u = 1 and 
T] = 1/4, at the critical point T = T^ = 2/ln(-\/2 + 1) (we set J = 1) and h = he = 0. 
Critical correlations do not develop a diverging length scale in the presence of the 
external space-dependent potential, i.e., at fixed v > 0. But a critical behavior develops 
in the limit w — )■ 0, described by TSS as in Eqs. ([2]) and (IHl). The trap exponent is given 
by 

9 = 16/31, (9) 

which can be derived by RG scaling arguments ^U\ taking into account the fact that the 
external potential couples to the order parameter. The static TSS has been numerically 
investigated in [10], here we focus on the dynamic TSS scenario of Eqs. dS]) and (jH]). 

We consider a purely relaxational dynamics (also known as model A of critical 
dynamics [12]), which can be realized by Metropolis updatings in Monte Carlo 
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simulations. The corresponding dynamic exponent z has been accurately determined in 
homogeneous systems by numerical equilibrium and off-equilibrium methods, obtaining 
z ~ 2.17.1j| According to the dynamic TSS scenario, in the presence of the trap we should 
not observe any diverging time scale r, but the critical dynamics is only recovered in 
the limit / — )■ oo, with a power-law diverging time scale ([6]). 

In our simulations we consider square lattices with —L < x,y < L, and trap 
potential V{x) = (x^ + y^)/l'^ (thus the trap is centered at the origin), while at the 
boundaries \x\, \y\ = L + 1 the spin variable is kept fixed: Sj = 1 (corresponding to 
Pi = 0). In the presence of the trap of size /, the lattice size L is chosen to have the 
spin variables effectively frozen at the boundary, making unnecessary the use of larger 
lattices to effectively obtain infinite volume results. See also below. The dynamics is 
provided by Metropolis updatings of the lattice variables using a checkerboard scheme 
(a time unit corresponds to a sweep of the lattice). We present results of MC simulations 
for several values of the trap size, up to / = 320. The typical statistics of our simulations 
is 10^ sweeps. 

We want to determine the power-law behavior of the time scale at T = Tc as a 
function of the trap size. For this purpose we estimate the autocorrelation times of the 
magnetization along the lattice axes, which can be easily computed. More precisely, we 
define 

[rlO] 

Mr= Yl ^iO + Soi (10) 

i=-lrlO\ 

where |_a:J = floor(a;) is the largest integer not greater than x, the distance r is measured 
in unit of the length scale l^ (to construct nonlocal observables compatible with the 
expected TSS behavior), and the subscripts of Sij indicate the coordinate of the lattice 
sites. In the following analysis we consider the values r = 1/2 and r = 1. We estimate 
the autocorrelation times of the correlators 

a(t2-tl) = (M,(ti)M,(t2))c (11) 

For this purpose we exploit the method of |35j which provides a substantial numerical 
advantage because it does not require extrapolated large-time estimates of the 
autocorrelation functions. We define 

'-<' + "^^>" ln|fiM/fi(t + „)] - <12' 

where n is a fixed integer number. A linear interpolation can be used to extend fr{t) to 
all real numbers. Then, for any positive x, we define an autocorrelation time r^,^ as the 

I Some of the most recent equilibrium estimates are z = 2.1667(5) from [5S], z = 2.168(5) from [5S], 
z = 2.1665(12) from [Uj, z = 2.172(6) from [2H]- Off-equilibrium results are z = 2.156(2) from ^, 
z = 2.153(2) from ^, z = 2.16(3) from ^, z = 2.1337(41) from [31], z = 2.160(5) from [33], 
z = 2.165(10) from [51]. The apparent small discrepancy between the two sets of results should not be 
significant. 
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Figure 1. The ratio ti/2(-^)/ti/2(oo) vs LI 



solution of the equation 



I r \ tXj I T 



(13) 



Tr^x are expected to show the same power-law behavior with increasing Z, i.e., t^^x ~ l'^- 
After some checks, we chose x = 2 as optimal value for our calculations, so in the 
following we present only results for a; = 2 and skip the subscript indicating the value 
of X. 

Let us first discuss the finite size effects due to a finite volume L^, which are expected 
to enter through the FSS variable Ll~^ [19]. This is confirmed by Fig. [H which shows 
the data of the ratio Ti/2(L)/ri/2(oo) between the autocorrelation times of M1/2 at a 
given L and in the L -^ 00 limit, as estimated by the result for the largest available 
lattice. The analysis of the FSS effects show that, within an accuracy of a few per mille, 
the data for Ll~^ > 5 can be effectively considered as L — )■ 00 data. 

In Fig. [2] we show the resulting infinite- volume estimates of the autocorrelation 
times of M1/2 and Mi, up to a trap size / = 320. They clearly show an asymptotic 
power-law behavior. The data for / > 48 turn out to fit the simple power law Tr = al^ 
with x^/d.o.f < 1, providing the estimate k = 1.120(3) (the results of the fits appear 
reasonably stable with increasing the minimum value of I allowed in the fit, while the 
data for / < 48 show small deviations from the simple power law). Analogous results 
are obtained by varying the value of x in the definition ( IT3l) . Dividing the exponent k 
hy 9 = 16/31, we obtain 

K/e = z = 2.170(6) (14) 

which is in good agreement with the best available estimates of z, for example [25] 
z = 2.1667(5). We have also performed fits to all data including scaling corrections, 
which are expected to be [H] 0(/~^), obtaining consistent results. 

§ This definition is based on the idea that, if the autocorrelation function C(t) were a pure 
exponential, i.e., C(t) = Co exp(— i/r), then T{t) = t for all t and thus t^ = t for any x. It 
provides a good autocorrelation time for any x, which converges to the integrated autocorrelation 
time Tint = ^ St^-00 C'(i)/C(0) for a; — > 00. See [35] for a through discussion. 
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Figure 2. Power-law behavior of the autocorrelation times T1/2 and ti. The lines are 
drawn to guide the eyes. 



In conclusion, our numerical results for the relaxational dynamics of the square- 
lattice gas model ([7]) show that the dynamic TSS scenario correctly describes the 
development of a critical dynamics in the presence of an inhomogeneous external 
potential which acts as a trap for particle systems. We expect that the validity of 
the dynamic TSS extends to other dynamic universality classes. 
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